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Using the monomer-dimer representation of strongly coupled U(N) lattice gauge theories with 
staggered fermions, we study finite temperature chiral phase transitions in (2 + 1) dimensions. A 
new cluster algorithm allows us to compute monomer-monomer and dimer-dimer correlations at 
zero monomer density (chiral limit) accurately on large lattices. This makes it possible to show 
convincingly, for the first time, that these models undergo a finite temperature phase transition 
which belongs to the Kosterlitz-Thouless universality class. We find that this universality class is 
unaffected even in the large N limit. This shows that the mean field analysis often used in this limit 
breaks down in the critical region. 
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Computing quantities in Lattice QCD with massless 
quarks is notoriously difficult. Most known algorithms 
break down in the chiral limit. For this reason questions 
related to the universality of chiral phase transitions are 
among the many questions that remain unanswered. It is 
often difficult to compute critical exponents sufficiently 
accurately to rule out all possibilities except one. 

The most useful simplification of lattice QCD occurs 
in the strong coupling limit which retains much of the 
underlying physics of QCD except for large lattice arti- 
facts. In this limit spontaneous chiral symmetry breaking 
and its restoration due to finite temperature effects have 
been studied using large N and large d expansions [1-3]. 
However, since these approaches are based on mean field 
analysis they cannot help in determining the universality 
of phase transitions. 

Interestingly lattice QCD with staggered fermions in- 
teracting through U(N) gauge fields can be mapped into 
a monomer-dimer system in the strong coupling limit [4] . 
These models contain an exact U(l) chiral symmetry, a 
remnant of the full chiral symmetry of QCD. When it 
was proposed, the monomer-dimer representation offered 
a new approach to study strongly coupled gauge theories 
close to the chiral limit from first principles. Unfortu- 
nately, this dream has remained unfulfilled until now. 
As in the weak coupling regime, most numerical simula- 
tions of the monomer-dimer systems have suffered from 
critical slowing down close to the chiral limit and hence 
have only allowed calculations with limited accuracy [5]. 

Recently, a cluster algorithm has been discovered to 
study these strongly coupled lattice gauge theories in the 
chiral limit [6]. This allows precision calculations in the 
chiral limit for the first time. As a first application of this 
new algorithm, in this article we study the finite temper- 
ature critical behavior in (2 + 1) dimensions. In agree- 
ment with expectations from universality, we find with 
very high precision the chiral phase transition to be in 
the same universality class as the Berezinski-Kosterlitz- 
Thoulcss (BKT) transition [7] . Our results are compara- 
ble to other known high precision spin-model studies of 
this universality class. We also show that the BKT tran- 
sition persists even in the large N limit, showing that the 



mean field analysis breaks down in the critical region. 

Although our results will primarily be of interest to 
lattice gauge theorists, it is possible that they may have 
a wider audience. We are basically studying monomer- 
dimer models and such models have a long history [8, 9]. 
They have been studied extensively in the context of sta- 
tistical and condensed matter physics. In the sixties, 
these models attracted a lot of attention [10, 11] when 
it was shown that the Ising model can be rewritten as 
a dimer model [12]. In the late eighties, these models 
became fashionable again in the quantum version [13] as 
a promising approach to the famous Resonating- Valence- 
Bond (RVB) liquid phase [14]. More recently the quan- 
tum dimer approach has again gained momentum since it 
was shown that the RVB liquid phase is actually realized 
on a triangular lattice [15]. This has raised hope that 
dimer models may also yield a theoretical understanding 
of the spin- liquid phase [16]. Inspite of these studies, as 
far as we know, a BKT transition has not been studied 
in the context of a dimer model. 

The partition function of the model we study in this 
article is given by ([4]) 



Z= [dUjdipdip exp(-S[U,ip,ip]) 



(1) 



where [dU] is the Haar measure over U (N) matrices and 
dip dtp specify Grassmann integration. At strong cou- 
plings, the Euclidean space action S[U, ip, tp] is given by 
the fermionic part 
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(ij) 



ipiU^tpj 



ipjUj^ipi 



m 22 ipiipi 



(2) 



where i refers to the lattice site on a periodic three dimen- 
sional cubic lattice of size L x ,L y along the two spatial 
directions and size L t along the euclidean time direction, 
(ij) represents the bond connecting the nearest neighbor 
sites i and j, Uuj\ are N x N unitary matrices associated 
with the bond (ij) and represent the gauge fields, ipi is 
an N component column vector and ipi is an N com- 
ponent row vector made up of Grassmann variables and 
represent staggered fermion fields at the site i. We will 
assume that the gauge links satisfy periodic boundary 
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conditions while the fermion fields satisfy either periodic 
or anti-periodic boundary conditions. The factors rj^ 
are the well known staggered fermion phase factors that 
depend on the coordinates z(or j). We will choose them 
to have the property that rj 2 ^ = 1 when (ij) is a spatial 
bond, and rfi^ = T when (ij) is a temporal bond. The 
exact form of T)/ij\ is not important at strong couplings 
since only rf^ appears in the final partition function. 
The real parameter T acts like a temperature. By work- 
ing on asymmetric lattices with L t « L and allowing 
T to vary continuously, one can study finite temperature 
phase transitions in strong coupling QCD [5]. 

The partition function given in eq.(l) can be rewritten 
as a monomer-dimer system given by 



^ = E IT>fe>) fc 

[n,b] (ij) 



(N-b {ij) )l 
b {ij) \N\ 



(3) 

and is discussed in detail in [4, 6]. Here rij = 0, 1, 2, N 



refers to the number of monomers on the site 



0. 1. 2, N represents the number of dimers on the bond 
(ij), m is the monomer weight, — 77?- arc the dimcr 
weights. Note that while spatial dimers carry a weight 

1, temporal dimers carry a weight T. The sum is over 
all monomer-dimer configurations [n, b] which are con- 
strained such that the sum of the number of monomers 
at each site and the dimers that touch the site is always 
N (the number of colors). The ' in the sum reminds us 
of this constraint. In this work we choose L x = L y = L. 

Comparing eqs. (3) and (1) we learn that zero 
monomer density corresponds to the chiral limit. The 
bipartite nature of the lattice can be used to distinguish 
every site j as either even or odd. If we define Oj = +1 
for an even site and aj = — 1 for an odd site, it is easy 
to show that when m = the action given in eq. (2) is 
invariant under U(l) chiral transformations, 



(4) 



The results of [6] show convincingly that this U(l) chi- 
ral symmetry breaks spontaneously at zero temperatures 
in two spatial dimensions. In this article we study the 
finite temperature critical behavior and show that it be- 
longs to the BKT universality class as expected. For this 
purpose we compute the chiral susceptibility Xs and the 
winding number susceptibility xw , both of which are of- 
ten computed in this context. The chiral susceptibility 
is equivalent to the spin susceptibility of classical XY 
models and is given by 



d 2 z 



VZ {dm) 2 



m=0 



(5) 



where V = L 2 L t is the lattice volume. It is easy to see 
that Xs is the integrated monomer-monomer correlation 



function. We define the winding number susceptibility 
by 



Xw = ^(W 2 + W 2 ). 



(6) 



where W x and W y are the winding numbers in the x 
and y directions respectively and are given by W x = 
(J2 z Vi b (ij))/ L x,&nd W y = {Y<i a i b (ii))/ L y In the def- 
inition of W x and W y the site j is chosen such that the 
bond (ij) is along the positive x and y directions respec- 
tively. 

In this study we fix L t and compute Xs and xw as 
a function of L. There are striking predictions for the 
large L behavior of these quantities if the phase transition 
belongs to the BKT universality class. If T c represents 
the critical temperature we expect 



and 



L 2- n (T) T<Tc 

u x < L 1 - 75 [log(L)] - 125 T = T C 
constant T > T r 



1/[2t7(T)] T < T c 

\u ; [2 + l/Iog(L/L )] T = T C 
T > T c 



(7) 



(8) 



in the large L limit. At T = T c the next to leading order 
correction is also shown since it can be important. The 
critical exponent r)(T) is expected to change continuously 
with temperature but remain in the range < f](T) < 
0.25 assuming the monomers carry a unit of U(l) charge. 
These predictions have been discussed in [7, 17-19] and 
have been used in [20, 21] earlier to demonstrate BKT 
behavior. In order to confirm these predictions we have 
computed Xs and xw using the algorithm discussed in 
[6], for lattice sizes ranging from L = 32 to L = 750 and 
the number of colors from N = 1 to TV = 32. We vary T 
with fixed L t — 4 to study the critical behavior. 

Let us first discuss our results for N = 1. In figures 1 
and 2 we plot Xs and xw as functions of L. We find that 
Xs fits extremely well to the form bL 2 ~ v when T < 1.0 
with 77 changing continuously in the expected region. In 
table I we show the fits along with their % 2 per degree 
of freedom (x 2 /d.o.f.). In the last column we give the 
results for l/(2xw) 011 the largest lattices that we could 
compute it. This value matches the value of r\ obtained 
by the fits. The figures and the fits also indicate that 
when T > 1.10 one is in the high temperature phase 
with a finite correlation length. 

When 1.02 < T < 1.06, we find that the values of 77 
do not match l/2%w even on the largest lattices while 
their values are close to a quarter as expected near the 
phase transition. There are noticeable finite size effects 
in xw ■ Assuming that these discrepancies are due to the 
logarithmic corrections we fit the data for Xs to the form 
bL 2 ~ n (log(L))~ 2r . In table II we show these new fits. We 
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FIG. 1: Plot of \s vs. L for various values of T. 
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FIG. 2: The plot of \w as a function of L for various values 
of T. The inset shows xVd.o.f vs. T for the fit xw = 2 + 
1/ log(L/Lo) which is valid only at T c . 



also fit the data for \w to the form (c+ 1 / log(L/Lo)) fol- 
lowing [20] . The last column of table II shows the values 
of l/2c obtained from the fit. We find that the loga- 
rithmic terms are unimportant for T < 1.0 as expected. 
However, they appear to affect the value of rj and l/(2c) 
(which is l/(2xw) in the thermodynamic limit) in the 
range 1.02 < T < 1.06. In particular the value of 2r is 
consistent with expectations at T — 1.04. Further, the 
new values of rj are closer to l/2c. Using the fact that 
Xw = (2 + 1/ log(L/L )) is exactly valid at T = T c we 
fit the data for \w to this form for various values of T. 
In the inset of figure 2 we plot x 2 /d.o.f as a function 
of T . Based on where the minimum occurs we estimate 
T c = 1.040(5). 

We have checked that \s and xw show similar evi- 
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0.2256(7) 
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0.2321(7) 
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0.2378(8) 


1.06 
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0.2421(5) 
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0.246(1) 
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0.413(1) 


0.2604(4) 
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TABLE I: Fits for x° = bL 2 ~ v . The last column gives the 
value of l/2xw on the largest lattices. 
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2r 


X 2 /d.o.f 


l/(2c) 


1.00 


0.398(5) 


0.222(5) 





0.2 
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0.249(5) 
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0.2583(5) 


1.10 


0.285(4) 


0.388(5) 


-0.12(2) 
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0.2831(5) 


1.14 
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0.569(6) 


-1.24(3) 
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TABLE II: Fits for X s = 6L 2 ""(log(L))- 2r . The last column 
gives the value of l/2c obtained from the fits to xw discussed 
in the text. 



dence for a BKT transition at higher values of N. Us- 
ing techniques similar to the JV = 1 we have com- 
puted T c for various values of N. We find that T c = 
0.708(6)7V + 1.40(4) - 1.07(4) /N fits our results very well 
for all values of N with a x 2 /d.o.f of 1.1. The detailed 
analysis of results at higher N and the dependence of 
non-universal critical properties on L t will be presented 
elsewhere. 

With regards to the N dependence of our results we 
find two interesting observations. When N approaches 
infinity it is often believed that a mean field analysis be- 
comes useful. Here we argue that there are exceptions 
to this lore. Consider the Mermin-Wagner-Coleman the- 
orem which shows that a continuous symmetry cannot 
break in two dimensions [22]. However, it can be shown 
that ?y = at N — ► 00 in the mean field approach, which 
means that the condensate is non-zero. Clearly, this can- 
not be true at any large but finite N. Witten has argued 
that when the symmetry is U(l) the large N analysis is 
still applicable since 77 ~ X/N at large N [23]. Our results 
for a fixed T = 1.0, shown in figure 3, do agree with this 
conjecture. Interestingly, as N becomes large and T/T c 
is held fixed instead of T, we find that rj ^ even in the 
large N limit. Figure 4 shows that 77 approaches an inter- 
esting function of T/T c as N becomes large. Extending 
this observation to QCD, we think that the t'Hooft limit 
(large N, small gauge coupling g, with g 2 N held fixed) 
may be quite similar [24]. 

In strictly two dimensions the monomer-monomer cor- 
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breaking of a global U(l) chiral symmetry in three di- 
mensions. The winding numbers W x and W y we have 
defined are exactly the magnetic flux in the x and the 
y directions obtained using the magnetic field defined in 
[26]. We interpret the divergence free magnetic field as 
the conserved current related to the E/(l) chiral symme- 
try. 
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FIG. 3: The plot of 77 as a function of N lor T — 1. The data 
fits very well to the form 77 = 0.169(6) /iV + 0.050(9)/jV 2 with 
X 2 /d-o.f.= 1.2. 
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FIG. 4: Plot of 77 as a function of T/T c . 



relations for N — 1 were calculated in [11, 25]. It was 
found that 77 ~ 0.5, which was recently confirmed again 
in [6]. This is not in the range < 77 < 0.25 we found 
above. We think a closer examination of this result in 
the context of BKT transitions would be interesting. 

Finally we note that the N = 1 model has been re- 
cently studied in [26] on various types of lattices, using 
another algorithm which is efficient in measuring dimcr- 
dimer correlations at zero monomer density [27]. The 
authors of [26] show the existence of a Coulomb phase 
in three dimensions on cubic lattices. Our results sup- 
port this observation indirectly, however we interpret the 
results differently. We use the connection of the dimer 
model with strongly coupled lattice QED to show that 
the long range correlations are a result of a spontaneous 
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